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Abstract 

In this paper we apply the formal Inverse Spectral Transform for 
integrable dispersionless PDFs arising from the commutation condi¬ 
tion of pairs of one-parameter families of vector fields, recently devel¬ 
oped by S. V. Manakov and one of the authors, to one distinguished 
class of equations, the so-called Dunajski hierarchy. We concentrate, 
for concreteness, i) on the system of PDFs characterizing a general 
anti-self-dual conformal structure in neutral signature, ii) on its first 
commuting flow, and iii) on some of their basic and novel reductions. 

We formally solve their Cauchy problem and we use it to construct 
the longtime behavior of solutions, showing, in particular, that unlike 
the case of soliton PDFs, different dispersionless PDFs belonging to 
the same hierarchy of commuting flows evolve in time in very different 
ways, exhibiting either a smooth dynamics or a gradient catastrophe 
at finite time. 


1 Introduction 

Waves propagating in weakly nonlinear and dispersive media are well de¬ 
scribed by integrable soliton equations, like the Korteweg - de Vries izu equa¬ 
tion and its integrable (2 -|- 1) dimensional generalization, the Kadomtsev - 
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Petviashvili [SU] equation. The Inverse Spectral Transform (1ST), introduced 
by Gardner, Green, Kruskal and Miura [2D], is the spectral method allowing 
one to solve the Gauchy problem for such PDEs, predicting that a localized 
disturbance evolves into a number of soliton pulses + radiation, and solitons 
arise as an exact balance between nonlinearity and dispersion H. 0.0.111. 
Soliton PDEs arise in hierarchies of commuting flows, and equations of the 
same hierarchy share the same multisoliton solution and similar behavior. It 
is known that, apart from exceptional cases, soliton PDEs do not generalize 
naturally to more than (2 + 1) dimensions; therefore, in the context of soliton 
equations, integrability is a property of low dimensional PDEs. 

There is another important class of integrable PDEs, the so-called disper¬ 
sionless PDEs (dPDEs), or PDEs of hydrodynamic type, including as distin¬ 
guished multidimensional examples the dispersionless Kadomtsev-Petviashvili 
(dKP) equation [6011671, describing weakly nonlinear and quasi one dimen¬ 
sional waves in Nature [HHDOlinT], the heavenly equation, arising from the 
the vacuum Einstein equations and the conformal anti-self-duality condi¬ 
tion for the signature metric in canonical Plebanski form [5D], the Duna- 
jski equation ra. an integrable generalization of the heavenly equation in 
which only the anti-self-duality condition in kept, the dispersionless 2D Toda 
(d2DT) equation [8l[T8l|66|, whose elliptic and hyperbolic versions are both 
relevant, describing, for instance, integrable "H-spaces (heavens) [H1[2T], in¬ 
tegrable Einstein - Weyl geometries [2D]- [2D], [DD], and playing a key role 
in the study of the ideal Hele-Shaw problem [DDllDDl IDTl IDT] : the Pavlov 
equation 1531, 071, 05, arising in the study of integrable hydrodynamic 
chains, and the Manakov-Sant ini system [12], giving a local description of any 
Lorentzian Einstein-Weyl geometry [T2], and including, as particular cases, 
the dKP and the Pavlov equations. dPDEs arise, more in general, in various 
problems of Mathematical Physics and are intensively studied in the recent 
literature (see, f.i., [Dl[7l[IlHI7l[221|28l[29l|3lH35l[Mll52Hp[^^ 

Since they arise from the condition of commutation [L, M] = 0 of pairs of 
one-parameter families of vector fields, implying the existence of common 
zero energy eigenfunctions: 

[L, M] = 0 ^ Lip = n't!: = 0, (1) 

they can be in an arbitrary number of dimensions [65]. In addition, due 
to the lack of dispersion, these multidimensional PDEs may or may not ex¬ 
hibit a gradient catastrophe at hnite time and, as we shall see in this paper, 
even in the same hierarchy there exist equations evolving in a smooth way 
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or into a gradient catastrophe at finite time. Their integrability gives a 
unique chance to study analytically such a mechanism and, also with this 
motivation, a novel 1ST for vector fields, significantly different from that of 
soliton PDEs |Il|2l|6l], has been recently constructed at a for¬ 

mal level, by Manakov and one of the authors of this paper (PMS), i) to 
solve their Cauchy problem ii) obtain the longtime behavior of 

solutions [HHIH], iii) costruct distinguished classes of exact implicit solu¬ 
tions [HHIT], see also |1], iv) establish if, due to the lack of dispersion, the 
nonlinearity of the PDE is “strong enough” to cause the gradient catastro¬ 
phe of localized multidimensional disturbances, and v) study analytically the 
breaking mechanisms [MIHSII^ . It is important to mention that the main 
difficulty to make the above 1ST for vector fields rigorous is associated with 
the proof of existence of analytic eigenfunction, motivated by the small field 
limit in the spectral parameter A. A proof of existence, but only for Im A > c, 
was found in |23] for the 1ST of the dKP equation; a complete proof exists 
so far only in the simplest case of the Pavlov equation, for which the whole 
1ST has been recently made rigorous |24j . 

In this paper we concentrate on the Dunajski hierarchy, whose elegant 
structure was enveiled in mu- More precisely, in §2 we consider some basic 
members of the Dunajski hierarchy, corresponding to the same Lax operator 
L: i) the system of PDEs |1] characterizing a general anti-self-dual conformal 
structure in neutral signature |T2] , ii) its first commuting flow, and iii) some 
of their basic and novel reductions. In §3 we construct the formal 1ST for the 
Dunajski hierarchy corresponding to the same Lax operator L, and we use 
it to solve the Cauchy problem for localized (in x, y, z) initial data. While 
the direct and inverse problems are the same for all the equations of this 
hierarchy, being associated with the same Lax operator L, the f-evolution of 
different equations of such hierarchy, ruled by different M operators, is con¬ 
siderably different, as we shall see in this paper, leading either to a smooth 
dynamics or to a gradient catastrophe of dKP type [44114^ . In §4 we con¬ 
struct the nonlinear Riemann-Hilbert (RH) dressing for the above equations, 
connecting also the RH data to the initial data of the dPDEs. In §5 we dis¬ 
cuss the constraints on the spectral data corresponding to the reductions of 
the first commuting flows of the hierarchy. In §6 we use the NRH problem to 
construct the longtime behavior of the solutions. In §7 we make some con¬ 
cluding remarks and we discuss interesting open problems to be investigated 
in future works. 

This paper is dedicated to the memory of S. V. Manakov. 
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2 The first members of the Dunajski hierar¬ 
chy and their basic reductions 

The Dunajski hierarchy is a basic example of hierarchy of integrable dPDEs, 
including the heavenly and the Manakov-Santini hierarchies as particular 
cases; see [ 1 ] for details and for an elegant characterization of it. 


2.1 The first two commuting flows 


The hrst and basic member of such hierarchy is the following system of three 
PDEs m 


Uxtl VyjUxy T UyllxX VxUyy fy, 

'^xti ^yz ^y')^xy T ^y'^xx '^x'^yy fxj 

fxt\ fyz ^y^fxy T ^yfxx ^xfyy 0; 

eqnivalent to the commntation [L, Mi] = 0 of the vector helds: 


L dz T ^dx T Ux fxd\i 

Ml = dt^ + \dy -\- Uy — fydx, 


( 2 ) 

(3) 


where 


u = udx + vdy. (4) 

It was recently shown [12] that there exist local coordinates {ti,z,x,y) snch 
that any anti-self-dual conformal structure in signature ( 2 , 2 ) is locally rep¬ 
resented by the metric 


g = dtidx + dzdy + Uydt\ — {ux + Vy)dtidz Vydz^ 


(5) 


where u and v satisfy eqnations ([ 2 |). 

Keeping L hxed and varying M, one obtains a hierarchy of commuting 
flows. The hrst commnting how reads 


'^xt2 if^z T 'dx')Oi Uyfx '^xyfd Uxx^ ) 

^xt2 (^2 T ^yfx "^xyfd Vxx^^y ( 6 ) 

fxt2 + Ux)^ fxfy fxy/d fxxOt-, 

where the helds a, /9 ,7 are dehned in terms of u, u, / through the equations 


Oix ^yz '^x^xy "^x^yy '^y'^xy “ 1 “ ’Uytlxx fy-) 

Px '^x'^xy '^x'^yy '^y'^xy “ 1 “ '^y’^xx ^fxi (7) 

^x ]yz “1“ '^xfy '^yfx' 
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Equations ([S]),® arise as the commutation condition [L, M 2 ] = 0 of L with 
the vector held 


M 2 — dt2 + {uy\ + a)dx + (A^ + Vy\ + (3)dy + {—fyX + 'j)d\, (8) 


whose coefficients increase by one their degree as polynomials of A, with 
respect to Mi. 

2.2 Some basic reductions 

The Dunajski hierarchy admits interesting reductions (see, hi. IHE]). The 
reductions of equations (E]) and (I6]),(I7]) we consider in this paper are of two 
types. Arising from the commutation conditions [L, M„] = 0 for one param¬ 
eter families of vector helds, the hrst and natural reduction considered here 
is the divergence free condition for such vector helds (the condition V • u = 0 
for the vector held X = dt + u ■ V), implying the constant volume condition 
for the associated dynamical systems. The second basic set of reductions 
explored in this work, giving rise to non-autonomous systems, are associ¬ 
ated with the condition that the so-called Orlov eigenfunctions of (L, M„) be 
polynomial in the spectral parameter A (see §5 for more details). 

Beginning with the system (EJ, the divergenceless constraint for L and 
Ml : 

'^X T Vy - 0, U - 6y, V - Ox (9) 

leads to the well known Dunajski equation [10] 



( 10 ) 


and to its vector helds Lax pair: 


L Qz T ^dx T Oxydx T Oxxdy fxdx, 

Ml = dt + Xdy + Oyydx “ 0 xyOy “ f yO X , 


( 11 ) 


reducing to the second heavenly equation 

Oxti ~ Oyz -\- OxxOyy — O^y = 


( 12 ) 


for / = 0. 
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The second set of reductions of system ([2]) lead to non-autonomous and, 
to the best of our knowledge, novel integrable nonlinear PDEs. 


The reduction 


leads to the system 


u + zf = 0, 


^xtl ^yz T T Vy^Vxy ^fyVxX VxVyy fx, 

fxti fyz T i^^fx T Vy^fxy ^fyfxx ^xfyy 0- 


The reduction 


v + kf = 0 , 


(13) 


(14) 


(15) 


plays a similar role and does not give anything new, leading to the system 


^xt\ ^yz (tlfy T Ux')^xy T fx^yy T 'U'yUxx fyi 

fxti fyz (flfy T Ux^fxy T tlfxfyy T Uyfxx 0, 

equivalent to fflTl) by the change of variables u v, z fi, x -H- y. 


(16) 


The combination of the reductions ([2]) and 0131) : 

Ux + Vy = u + zf = 0 u = z9y, V =-z9x, f =-9y, (17) 

leads to the following scalar non authonomous generalization of the second 
heavenly equation; 

9xti — 9yz + z9xx9yy — z9‘^y = 0, (18) 

while the combination of the reductions ffT3|) and flT^ 

u + zf = V + tif = 0, (19) 

leads to the scalar dPDE 


fxt fyz T zfxfxy T tifxfyy zfyfxx tlfyfxy 0. 


( 20 ) 


Concentrating now on equations mM, the divergenceless constraint for 
the vector helds L, M 2 , expressed now by the equations 

Ux + Vy = ax + {/3-f)y = 0 u = 9y, V =-9x, a = py, 13 = f - px {21) 
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leads to the system 


((^Xt2 ~ {dz + 0xydx ~ dxxdy)p + (^xyf^ + (^z + (^xydx — 0xxdy)f 
fxt 2 fxfy “1“ (^z "h ^xydx ^^xx^y')'^ "h {^f Px^fxy "h Pyfxx 0) 


where 


px 

lx 


+ ^yz 


+ ^ly 


^xx^yy f 0, 


‘^^xyfxy ^xxfyy dyyfxxi 


reducing to the second member of the heavenly hierarchy 


0, 


(22) 

(23) 


^xt2 if^Z “1“ ^xydx ^xx9y') Py 

Px "t~ dyz “1“ ^xy ^xx^yy 0 

for / = 0 and, consequently, 7 = 0 . 

The reduction (IT^ leads instead to the system 


(24) 


'^xt2 (^z zfxdx T Vxdy^f^ fx'l^y ^'yVxx f^'^xyy 

fxt 2 T fxfy T (^dz ^fxdx T Vxdy^'y T f^fxy T fxx 0, 

where a = and 


(25) 


px {,9z ^fxdx '^x^y)'^y “1“ ( ^fy^x 1~ '^ydy^'^x “t“ ^fx") 

lx (^ 2 : ^fx^x 1~ "^X^y) fy ( ^fyQx 1~ "^ydy^ fx'j 


(26) 


and the combination of the above two reductions leads to the equations 


^xt 2 dy^xy Pz T zOxyPx zOxxPyy 

Px T ^yz T zO^y zOxx^^yy 0, 

where flTT)) hold and p = z~^p^ a = zpy^ (3 = —9y — zpx, 7 = Py. 
At last, the (less obvious) differential reduction 


(27) 


(1 - 2t2fy)v^ = 2 ^ 2 (/^ + 'ttrfx): 


(28) 


leads to the system 


(1 2 /^ 2 /*^) ^yfx I' '^xyP I' "^xx^ I' 1 ^x^x^ ^^2{,fz I' "^xfx^^y) 

(1 ^t 2 fy^(^fxt 2 "h fxfy H“ fxyP H“ fxx^ H“ Iz d* ^xlx^ ^t2(^fz H“ ^xfx')ly) 

(29) 
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where 


(1 2t2/'y)(c^a; “t“ Uyz UxUxy UyUxx fy) 2^2 [(O^ “1“ Uyfx)Uxy i^fz "h Uxfx^Uyy^, 

(1 — 2t2fy)(^ = 2 / + 2t2{ft2 + Q^/x), 

(1 2 ^ 2 / 2 /) (O'x fyz '^xfxy '^yfxx') 2^2 [(/z Uxfx')fyy “1“ '^yfx^fxy]- 

(30) 

The combination of reductions ffT3|) and fl28l) leads to the dPDE 

(1 - 2(2/„)(/j,b + /i/„ + fs,yl3 + 7, + - 2/i7i) = 2 * 2 ( 2 /^ - fy)jy, 

(31) 


where a = z'y and 

(1 — 2t2fy)/3 = 2/ + 2t2{ft2 + 

(1 - ‘it2fy){Zx - fyz + !:fxfxy “ ^fySxx) = 2i2|(/, - zf^)fyy - (t - zfxfy)fxy]. 

(32) 


At last, the combination of reductions fl2Tll and fl28ll leads to the PDEs 

2^2 [^Xt 2 ~ (^2 + ^xydx — dxxdy)p + ^xyf] — = 0, 

(1 — 2t2fy)0xx + ‘2t2{fz + (^xyfx) = 0, (33) 

Px + Oyz + — OxxOyy + / = 0 , 

where dH]) holds, and 7 = {6xyf)y - {0yyf)x - 0 xyl{ 2 t 2 ). 

The solution of the Cauchy problem of all the above reductions, in terms 
of a nonlinear RH problem, will be presented in §5. 

3 1ST for the Dunajski hierarchy 

In this section we apply the 1ST method introduced in [lOHll] to construct 
the formal solution of the Cauchy problem for the systems ([2]) and 0-0 
in (3 + 1) dimensions and for their reductions, within the class of rapidly 
decreasing real potentials u, u, /: 

M, U, / —)■ 0, + y‘^ + z^) +CX0, 

u,u,/gM, (x, y, ;s) G M^, fj > 0; (34) 

here tj, j = 1,2 are interpreted as time variables and x, y, z as space variables. 


3.1 Basic eigenfunctions 

A basic role in this 1ST theory is played by the real dost eigenfunctions 
V±{.x,y,z]\) G for A G M, the solutions of L0± = 0 dehned by the 
asymptotic: 

/ ^±o{x,y,z]X) \ ^ / A \ 

ip±{x,y,z]\) = \ ip±i{x,y,z]X) -)-^= ^ = x - Xz , ^±cx). (35) 

\ip± 2 {x,y,z-,X) J \ y / 

These dost eigenfunctions are intimately connected to the dynamical system 
(in the time variable z): 

/ -fx{x,y,z) \ 

— = X + u^{x,y,z) , (36) 

^ V v^{x,y,z) ) 

for the unknown x{z) = {X{z),x{z),y{z)Y in the following way. 

Assuming that the potentials u, v, f be smooth and sufficiently localized 
functions of x, y, z, it follows from ODE theory that the solution 

/ A{z;xo,yo,zo,Xo) \ 

x{z) = \ X{z]Xo,yo,zo,Xo) (37) 

\ Y{z-,Xo,yo,Zo,Xo) / 

of fl36|) . satisfying the initial condition {X{zQ),x{zQ),y{zQ))'^ = (Aq,xq,|/ o)^, 
exists unique, and it is globally dehned for real with the asymptotic states 
>^±,x±,y±: 


A ~ A±(a;o,2/o,2:o, Ao), 2 ;-)■ ±oo, 

X ~ zX±{x(i,yo,Zo,Xo) + x±{xo,yo,Zo,X(i), z± 00 , 
y ~ y±{xo,yo,zo,Xo), z± 00 . (38) 

These asymptotic states X±{xo,yo, Zq, Xq), x±{xo,yo, Zq, Xo),y±{xo,yo, Zq, Xq) 
are constants of motion for the dynamical system fl5B]l when the point (xq, yo, 
Zo, Ao) moves along the trajectories; therefore they are solutions of the vector 
held equation 

/ A±(x,2/,z,A) a 

L I x±{x,y,z,X) = 0 (39) 

\ y±{x,y,z,X) ) 
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and, due to fl5S]) and fISSD . they coincide with the real dost eigenfunctions 

^±o{x,y,z]\) = \±{x,y,z,\), 

ip±i{x,y,z]\) = x±{x,y,z,X), (40) 

^± 2 {x,y,z]\) = y±{x,y,z,X). 

A crucial role in the 1ST for the vector held L is also played by the analytic 
eigenfunctions y, z; A), the solutions of L'lp^ = 0 analytic, respectively, 

in the upper and lower halves of complex A plane, satisfying the asymptotics 

'ip^{x,y,z]\) ^ {x^+ y^ + z^) ^+ 00 . (41) 


They can be characterized by the following integral equations 

^l)'^{x,y, z; X) + / dx'dy'dz'G^{x — x',y — y', z — z'] X)' 

[u^/{x\y',z')d^' + v^[x\y\z')dy> - f^‘{x',y', z')dx]i^'^{x',y'X] A) 

= f, (42) 


for the analytic Green’s functions of the undressed operator {d^ + A^a,): 

^(y) 


G^{x,y,z;X) = ±— 


27ri[a; — (A + ie)z\' 


(43) 


Consequently, for |A| S> 1: 


where 


r±( AA c-*, Q{x,y,z) 

^ (x, y, A) = ^ +--+ 0(A ), 


f{x,y,z) 

Qix, y,z)= I -u{x, y, z) - zf{x, y, z) 

-v{x,y,z) 


(44) 


(45) 


We observe that the analytic Green’s functions (1431) exhibit the following 
asymptotics for z ±cx) 


G^{x — x\y — y\ z — z']X) —)■ ± 


Ky) 


2vri[^ - =F ie] ’ 


z —)■ +CX), (46a) 
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G^ix — x',y — y', z — z'-^X) ±—77-tt; z ^—00, ( 46 b) 

^ ’ 27ri[^ - ± k] ’ ’ ^ ^ 


where ^ = x — Xz, ^' = x' — Xz', entailing that the 2; —>• +cxd asymptotics of 
-0^ are analytic, respectively, in the lower and upper halves of the complex 
.^-plane, while the —)■ —00 asymptotics of -0^ are analytic, respectively, in 
the upper and lower halves of the complex ,^-plane. 

3.2 Scattering and Spectral data 

The z (time) scattering problem for the ODE system fl36|) allows one to 
construct the scattering data a(.^), dehned by the —)■ +cxo limit of 

lim 0_{x,y,z]X) = S{i) =i+a{i). (47) 

z^-\-oo 

The hrst part of the direct problem is the mapping from the real potentials 
/, u and V , functions of the three real variables (x, y, z), to the real scattering 
vector a dehned in fl47)) . function of the real variables ^ = {X,^,y) . 

Together with the above scattering data, one dehnes also spectral data 
in the following way. The ring property of the space of eigenfunctions allows 
one to express the analytic eigenfunctions in terms of the dost eigenfunctions, 
used as a basis for such a space, if A G M: 

ij+{x,y,z-,X) = lCt{0-{x,y,z-,X)) = IC-{(f+{x,y,z-,X)), (48) 

ij-{x,y,z;X) = JCZ{(p-{x,y,z-,X)) = JCX{ip+{x,y, z] X)), (49) 

where 

/C,“(0:=e + ;e^(0, a,b = ± (50) 

The 2 ; —)■ —00 and 2 ; —>• +cxo limits of respectively the hrst and second 
equalities in fl48|l and fl49|l imply 

lim =^+ (0, lim = ^ + y^(f), (51) 

z^—oo z^-\-oo 

Therefore the analyticity properties of ^|J^ established above imply that ;;!^(0 
and A^(0 analytic in Im ^ > 0 decaying at ^ cx) like 0(^“^), while 
Xl(.^) and xA(0 analytic in Jm <0 decaying at cxd like 0(.^“^). 
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In addition, taking the z —)■ +oo limit of the second of 048 p . one obtains the 
following equation: 

(f + - x;(0 =0, A e M, (52) 

that must be viewed as three “linear scalar Riemann-Hilbert (RH) problems 
in the variable with the given shift cf(0” for the unknowns X- ci.nd x+ 
(see, f. i., [IH] for the associated theory). Such a RH problems with a shift 
are equivalent to the following three linear Fredholm equations [19] 

x;.(0 -7^ [ K{i,e)x+M')d^' + = 0, j = 0,1,2, A e R, (53) 

^ 27r* Jr ^ 


where 


K(f = ^■S'(g')/ag' _ 1 

s(i'}-s(^ «'- 5 ’ 


MM) = -MM) + MPM. /r §^fM')di'. J = 0,1, 2, 

e = (A,e,2/), e' = (A,e',2/)- 

We remark that, if m, n, / G R, then the vector helds are real for A G R, and 


0± G R^, "0” = M^ A G R, 

a 6 K®, C- = CJ, x" = xi 


(64) 


Summarizing, the direct problem consists of the following steps: i) given 
the initial data {u{x, y, z, 0), f (a:, y, z, 0), f{x, y, z, 0)), one constructs the scat¬ 
tering data a(^, 0) from the solution of the ODE system (l36jl : ii) known 
(j(.^, 0), one solves the RH problem with a shift fl5^ constructing the spec¬ 
tral data M and M- 


3.3 Inverse problems 

First inversion. Known the spectral datum ;vt(^), an inverse problem can 
be constructed from the hrst of equations (HHp . observing that 

(55) 
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where 


h*s(A) - 


(56) 


S(A') 


A' — (A ± ie) 


dX' 


are the (±) analyticity projectors to the upper and lower halves of the com¬ 
plex A plane. Since P^/(A) = ±^Hxf{X) + Xf(^X)^ equation (l55|) is equivalent 
to the inverse problem formula 

(p_(x, y, z; A) + y, z; A)) + y, z; A')) = (57) 


where and 3";)^ are the real and imaginary parts of i.e., = 

, and Hxg{X) is the Hilbert transform operator 

Hxg{X) = -P.V. [ (58) 

^ Jr A — A 

Equation (157)) must be viewed as a vector nonlinear integral equation for the 
unknown dost eigenfunctions (p-. 

Once (p- is reconstructed from X- solving the nonlinear integral equation 
(1571) . f,u,v are hnally reconstructed from (see [M] ) 

/ = -^/R^xii(<^-(2:,2/,2:,A))dA, 

^ (59) 

u = -zf + y 4 Sxl2((^-()r, y, z, X))dX. 

To obtain the reconstruction formulae (EH]), one subtracts equations fHHjl and 
(jUj), and then uses the reality constraints flSTl) : 


- -0 = ^{p_) - X-{0-) = 2iQx-{0-), 
together with the analyticity properties of -0^, to get 


Equations 


= ^+- , V -n • ^ 

TT Jr A' - (A ± te) 

follow then from fHSjl . 


-dX'. 


( 60 ) 


( 61 ) 


Second inversion: the nonlinear Riemann-Hilbert problem. An al¬ 
ternative inversion is based on the nonlinear Riemann-Hilbert (NRH) inverse 
problem on the real A axis: 

^+{X) = n (i^-{X)^ , A e M (62) 
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or, in component form: 

/ Vi(A) \ / Ko (<^„-(A),t^r(A).Vf(A)) \ 

V>t(A) = K, (4-(A),^,-(A),V>2-(A)) , A€R, (63) 

V V>2*(A) / V K2 (,*„-(A),,*r(A),*(A)) / 

where the RH data 7^(C) are constrncted from the spectral data Xb(0 via 
algebraic manipnlations, eliminating 0^ from the hrst of equations fHSjl and 
The solutions (A), V’f (A),t/’^(A))^ G are analytic re¬ 

spectively in the upper and lower halves of the complex A plane, with the 
following normalization in the neighborhood of A = cxo: 

-?/>^(A) = A -h 0(A“^), V’^(A) = --^A + X + 0(A“^), 

V’^(A) = 2 / + 0(A"^), |A| > 1 , (64) 

and the potentials M,n, /, are reconstructed from (1441) : 

/ = lim A (V'^(A) - A), 
u = —zf — lim A ('0i^(A) — X + 2 ;A), 

V = hm X (ijjf (X) — y) . (65) 

3.4 Evolution of the scattering data 

In all formulas of the direct and inverse problems we have deliberately omit¬ 
ted the time variable, appearing just as a parameter. Now it is the time to 
introduce it, and, as we shall see, i) the evolution of the data is given by ex¬ 
plicit formulae (and this is the main justihcation for the introduction of the 
1ST); ii) such formulae imply a substantially different evolution for different 
members of the hierarchy; iii) the ti and ^2 flows commute on the level of the 
data, and this is the simplest way to prove that the corresponding flows () 2 |) 
and 0.0 commute as well. 

It turns out that, as f,u,v evolve in time according to ([2]), and, respec¬ 
tively, 0 . 0 . the time dependence of the three components of the data S, 
ICl and TZ are defined by the same PDEs: 


+ X^d^^)Dj = 0 , j = 0 , 1 , 

(dt. + A"%)Z^2 = D^, 
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whose explicit solutions read 

^.(Co, Cl, C 2 , tn) = D.iCo, Cl, C 2 - tnCo , 0), J= 0 , 1 , 

-D2(Co, Cl, C 2 , tn) = -D2(Co, Cl, C 2 — ^nCo , 0) + ^n(-Do(Co, Cl, C 2 “ ^nCo , 0))"', 

(67) 

where n = 1 and, respectively, n = 2. For an arbitrary n G N’*' in fl 66 |) and 
dn?]), one generates the time dependence of the spectral data of a sequence 
of commuting flows of the Dunajski hierarchy. 

Equations fl 66 p follow from the observation that one constructs, from 
the dost Lp- and analytic "0^ eigenfunctions of L, the common dost (/>_ and 
analytic eigenfunctions of both L and Mn operators via the formulae: 

(j)i{x,y,z,tn]X) = (p-o{x,y,z,tn;X), 

(p2ix,y,z,tn;X) = ip-i{x,y,z,tn;X), 

(j) 3 {x,y,z,tn] X) = ip- 2 {x,y,z,tn;X)-tn{(p-o{x,y,z,tn;X))'^ (68) 

'^Q{x,y,z,tn-,X) = ijQ{x,y,z,tn,X), 

'^^{x,y,z,tn-,X) = ijf{x,y,z,tn,X), 

^f{x,y,z,tn;X) = ^/J^{x,y,z,tn;X)-tnii/Joix^y^z^tn^X))"^, (69) 

using arguments developed, hi., in [ 12 ]. 

It is straightforward to verify that the flows fl 66 l) for different n commute. 


4 The nonlinear RH dressing 

Since the RH data 77(C, tn), n = 1 , 2 satisfy equations fl 66 l) . fl^ and the ana¬ 
lytic eigenfunctions of L are connected to the analytic eigenfunctions 
of L,Mn, n = 1,2 through the equations dSHD, it follows that the NRH for¬ 
mulation of the inverse problem, corresponding to the analytic eigenfunctions 
of both operators L, Mn, n = 1,2 reads as follows. 

Consider the NRH problem 

^+{X) = n(^^-{X)y AeR ( 70 ) 

or, in component form: 

/ 'I'+(A) \ / Ko ('I'„-(A),'I'r(A).l'2-(A)) \ 

®+(A) U Ki(®„-(A).®r(A),'I'2-(A)) . AeK, (71) 

V ('l'o~(A).«r(A),'I'2-(A)) J 
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where TZj{C) = iZj{(,0), the solutions ^^(A) = (\I/q (A), (A))^ G 

are analytic rispectively in the upper and lower halves of the complex A 
plane, with the following normalizations in a neighborough of A = oo: 

vI/f(A) = i/,(A) + 0(A-'), j = 0,l,2, (72) 

where, respectively, 

z/o(A) = A, iyi{X)=x-Xz, iy 2 {X) = y - Xh, (73) 

and 

I/O (A) = A, z/i(A)=x-Az, i/ 2(A) =-A^t2 + 1/- 2/^2- (74) 

Then 4/^(A) are vector eigenfunctions, respectively, of the vector helds L, Adi: 

= 0, and of the vector helds L,M 2 '- = 0; 

furthemore, the helds /, ti, n, reconstructed by, respectively, the following 
formulae: 

/ = lim A (T^(A) - A), 

A —y oo 

u = —zf — lim A (Ti‘^(A) — x + ; 2 A) , (Y 5 ) 

V = -tif - lim A (T^(A) - y + hX) , 

A^oo 

and 

/ = lim A (A) - A), 

A^oo 

u = —zf—\imX{'Alf{X)—x + zX), 

= (1 - 2 t 2 fy)~^ + Uxfx) - 

^3, lim A(T^(A) + A^f 2 -|/ + 2/f2)), (76) 

are solutions respectively of the nonlinear systems of PDEs ([2]) and 0 , 0 ^ 
In addition, the hrst few terms of the A large expansions of the above 
eigenfunctions of (L,Mi) and (L,M 2 ) read, respectively 

T±(A) = A + /A-i + A-2 + 0(A-3), 

T±(A) = -zX + x- {u + zf)X-^ + (5f^A-2 + 0(A-3), (77) 

+ 2/~ ^ + <52^A ^ + 0(A ^), 
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where 



+ ^xdx + Vxdy)f, 

\ + uA + v^dy){u + zf), 

^?x = + '^xdx + Vx:dy){v + fi/). 

(78) 

and 

T±(A) = A +/A-1 + A-2 + 0(A-3), 

^±(A) = -zX + x-{u + zf)X-^ + + 0(A-3), 

^±(A) = -t 2 A 2 + y- 2 ft 2 + + 0(A-2), 

(79) 

where 

(*^2 )a: 2^2(y^ T ^xfx T ^xfy) W- 

(80) 

The coefficients 52^^ and 52^^ satisfy also the compatible equations 



= (^ii + '<^ydx + Vydy){v + tif), 

{S 2 ^)y = 2 ^ 2(7 + Uyfx + Vyfy) - Vy. 

(81) 


The proof of these results is standard, in the NRH dressing philosophy 
(see, f. i., [H] for details). The NRH dressing for the Dunajski hierarchy was 
presented in [4], but no connection to the initial data was given. 

We remark that the NRH problem flTOll is characterized by the following 
nonlinear integral equations in the spectral variable A: 

^ L ^ 

where Rj{() = 7lj{() — Q, j = 0,1, 2, and the reconstruction formulae for / 
and u in terms of the RH data read 

while that for v is respectively 

»’ = -*i/+2l7/i.«=(®VA))<iA, (84) 

and 

V^ = {1- 2t2fy)~^ (2t2{fz + u^fx) + , (85) 
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for equations 02]) and (E]),©- 

The reality constraint for the RH data is a nonlinear formula: 

n = I ve e cl (86) 

It is possible to show that, if IZ satishes fl86l) . then 

f-(A) = T+(A), AgM (87) 

and -UjU, / G M (see, f.i., [H] for the proof). 

We end our remarks noticing that, in the normalization fl71]) correspond¬ 
ing to the second equation of the Dunajski hierarchy, the dependence on the 
variable y is through the combination y — 2/^2- Therefore the eigenfunctions 
T” and, through the reconstruction formulae (I83|) . (l8^ . also the solutions 
f,u,v, depend on y through the combination y — 2/^2- This is conceptually 
similar to the case of dKP |13], and provides the spectral mechanism for 
breaking of smooth and localized solutions of equation ([6]),([7]) at finite time 
^ 2 - Such a mechanism is absent in equation ([2]), for which a smooth evolution 
is expected. 

5 Constrained spectral data and reductions 

In this section we show how the eigenfunctions and the spectral data are 
constrained for the reductions of §2.2. We have already investigated in §3 and 
§4 the reality reduction. We concentrate now on the divergenceless constraint 
and on the constraint that the Orlov eigenfunctions are polynomial in A. 


5.1 The divergenceless constraint 

Consider a basis of three independent eigenfunctions V’j, 3 = 0,1, 2 of L and 
construct the 3x3 Jacobian matrix 

M = () (88) 

V d[XQ,Xi,X2) J 

(of components Mij = dtjji/dxj, i,j = 0,1,2), and its determinant J('0) = 
detM, where Xq = X, Xi = x, X 2 = y. Then the following identity holds 
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true: 


LJ 


[u. 


'^xy) Ji 


(89) 


implying that, if L is divergenceless, i.e., if satisfy the constraint ([9]), 
then J('0) is also an eigenfunction of L: = 0. 

In the direct problem philosophy, it is possible to show that, if u, v satisfy 
the divergenceless constraint ([9]), then the spectral data S,]C,7l satisfy the 
following constraint 


det 


fd{Do^D^M\ 

\ ^(Co,Cl5C2) / 


1 . 


(90) 


To show it, consider as a basis the dost solutions 0_. Then LJ(0_) = 
0 and J(0-) ^ 1 as z ^ —oo. Since 1 is an exact eigenfunction, then 
J(0_) = 1. Analogously, ~ 1 as A ~ oo; therefore are analytic 

eigenfunctions going like 1 at A ~ oo; then = 1. Evaluating J{4>-) = 1 

at 2 ; = oo, we infer the constraint fl90|) for the scattering data S. In addition, 
from equations flTSl) and flT9)) . it follows that 


J(T±) 


det 


I 9(Co,Ci.C2) J 


(91) 


implying the constraint fl90|) for the spectral data /C. At last, from the NRH 
problem fl7T]l . it follows that 


J(^+) = det 


V ^(Co,Ci)C2) / 




(92) 


implying the constraint fl90|) for the RH data TZ. 

Viceversa, in the inverse problem, if the RH data satisfy the constraint 
fl90|) . i.e., the NRH problem flTOjl is a volume preserving mapping, then m,u, 
reconstructed via (I75il . (|76ll . satisfy the divergenceless constraint (|9]). The 
proof is standard (see, f.i., |1]): if the RH data 7^(C) satisfy the volume pre¬ 
serving constraint ([90]), then, from (|92|), it follows that J(T+) = J(\['~), pro¬ 
viding the analytic continuation one of the other. Therefore J(tl/+) is entire in 
A, and since T(T'*') ^ 1 as A —)■ 00 , it follows that that T(T’^) = J(T“) = 1. 
As a consequence of it, the coefficient —{ux + Vy) of the 0(A“^) expansion of 
J(T“*') must be zero, and the constraint ([9|) follows. □ 
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5.2 Polynomial Orlov eigenfunctions 

Let g{X) be an eigenfunction of L characterized by the formal expansion 

N 

g{X)^ 9nX^ |A|»1, iV>0. (93) 

n=—oo 

Then the coefficients gn satisfy the recursion relations 
0 = gN,x, 

0 gN—i,x T T Ux^gN) 

0 gn—i,x T l)/'x 5 ^n+l T (^dz T Ux)gni ^ ^ XJ. (9"^) 

This three term recursion implies that, if g-i = 0, then gn = ^ for any 
n < —1. Consequently g-i = 0 is an admissible constraint for the dPDEs 
associated with L, and corresponds to the existence of polynomial eigenfunc¬ 
tions g{X). 

Let us apply this mechanism to the Orlov-type eigenfunctions and 
of the first two members of the hierarchy. The reduction ffT5]l is the g^i = 0 
condition for the eigenfunctions (see d77|) . d79|) h therefore it corresponds 
to the following elementary constraints on the eigenfunctions and data: 

= X — Xz Ri = 0 (95) 

for both dynamics (E]) and (I6]),([7]). The reduction ffT5|) is the g-i = 0 condition 
for the eigenfunctions Tf for the dynamics (EJ, and it corresponds to 

= y — Xti i ?2 = 0. (96) 

The g^i = 0 condition for the eigenfunctions for the dynamics 0.0 
corresponds instead to the differential reduction (E5]) (see fl79]) . (l5U]) ): it cor¬ 
responds to the constraints 

= —X^t 2 y — 2/^2 R 2 — 0 - (9”^) 

Using the above constraints on the RH data, we can characterize the solu¬ 
tions of the reductions presented in §2 in terms of NRH problems. 
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1. Consider the NRH problem and let the RH data satisfy the volume 
preserving condition 


det 


a(7^o(C),7^l(C),7^2(C)) 


= 1 ; 


^(Co, Cl) C 2 ) 

then f and 6, defined by 

= -^/ + 2 k Im K"(A). 'I'r(A). (A))dA, 


(98) 


(99) 


solve the Dunajski equation ma and equations (dl) for, respectively, the 
normalizations ((7^ and 

If, in addition, Rq = 0, then f = 0, d'g = A, |y7| ) becomes the 2-vector NRH 
problem 

vl/+(A) = 7^l (a, vI>-(A), d/ 2 -(A)), d/+(A) = 7^2 (a, vI/r(A), vI/-(A)), (100) 

the RH data satisfy the constraint 

a(7^l(C),7^2(C)) 


det 


5(Ci,C2) 


= 1 , 


and 


«.= 2 k 1.1*1 (^.«'r(A).'I'?(A))rfA 


( 101 ) 


( 102 ) 


solves the heavenly equation / flTI) and equations [2^ for, respectively, 
the normalizations (fTT]) and (f?^. 


2. Let Ri = 0; then u = —zf, \I/^ = x — Xz, |y7| ) becomes the 2-vector 
NRH problem 

d/+(A) = noho{X),x- Xz, '^f{X)\ 
d/2+(A) = 7^2(v[/o (A),a; - Az, ^/^-(A)), 
and f, V, constructed via 


(103) 


f = -^i W)dX, 

V = -tj + ^J^R 2 (^%iX),x - Xz, d'^(A))dA. 


(104) 
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and via 


f = 


‘2'ni 


Roi'^o 


= {l-2t2fy) ^ (2t2{fz- zf^) + j R2 {^q{\),x - \z,^ 2 {X)^d\ 


(105) 


for, respectively, the normalizations ((7^ and 


solve equations o and 

0 - 

If, in addition, the RH data satisfy the volume constraint 


det 


a(^o(a^)) 

^(Co, C2) 


= 1, 


(106) 


then [Tl\) are satisfied and 


6y = - 

271 i 


Roi^o {X),x - Xz ,^2 {X))dX, 


(107) 


solves equations / fl^) and |77| ) for, respectively, the normalizations (fyg| ) and 

O- 

3. Let R 2 = 0; then v = —tif, = jy 2 ^\X) for the normalization [7^ , 
and V satisfies the constraint (d^j and d'f = r'^\X) for the normalization 
(f^j), where of course 


= y -Xti, ^'f^(A) = -XH 2 + y- 2 ft 2 . 
The NRH problem ( [y7| ) becomes the 2-vector NRH problems 

>K+(A) = Ki(®,-(A), ®r(A). 4"’(A) ).n = 1, 2 


(108) 


(109) 


for the normalizations and and 


/ = -ife /, (®o (A). 'i'r(A). •^Pix))dx, 

= -^/ + ik /. ft («'o'(^). «'r(A). pir\x})d\ 


( 110 ) 
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are solutions of /[lU) and for respectively n = 
If Ri = R 2 = 0, then we have also u = —zf, 
becomes the scalar NRH problems 


1 and n = 2. 

= X — \f, and d77|j 


= 7^o(^o - A^, n = l,2 (111) 


for, respectively, the normalizations and At last, 

f = - Xz,4'^\X)^dX, n = l,2 

solve respectively the scalar dPDE 1^) and equation ^3l\). 

If R 2 = 0 and the RH data satisfy the volume constraint 


det 


d{Co, Cl) 


= 1 , 


( 112 ) 


(113) 


then f and 9, defined by 

/ =-jfe/. «o(l'„-(A), l-rlA). 4“’(A))<iA, 


(114) 


solve equations 


6 The longtime behaviour of the solutions 

As in the classical 1ST method, the two inverse problems for dPDEs are 
expressed in terms of integral equations in the spectral variable A, and the 
dependent and independent variables of the nonlinear dPDE appear there 
just as parameters. Therefore they are convenient tools to construct the 
longtime behavior of the solutions of the dPDEs, and this section is devoted 
to this goal, concentrating on the NRH problem. 

We start with the nonlinear integral equations (|82il . characterizing the 
NRH problems associated respectively with the dynamics ([2]) and (I6]),(I7]), cor¬ 
responding respectively to the normalizations (|73il and (Clj), and we rewrite 
it in a more convenient form 

= . = 0.1.2 (115) 
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in terms of the functions defined by 

d),(A):=v[/T(A)-i/,(A), J = 0,l,2. (116) 

Consider first the case (175]) . and 

fi > 1, x = viti, y = v2ti + y, z = V3ti, y,Vi,V2,Vs = 0(1) (117) 


Then equations 01151) read 

^ (a' + 1>o(A'). -V,(X - + 1>.(A'). 

-(A' - !>2)*i + S+42(A')), i = 0,1,2 


For t >> 1, the decay of the integral is partially contrasted if V 2 = or, 
equivalently, on the surface 

y--ti = y (119) 

z 

of the {x, y, z) space, on which the second and third argument of Rj (j = 
0,1, 2) in OllSp have the same linear growth in ti. On such a surface, since 
the main contribution to the integrals occurs when A' ~ U 2 , h is convenient 
to make the change of variables y! = (A' — V 2 )ti, obtaining, for ti » 1: 


—/i'+ ^ + $ 2(^2 + ^)), j = 0 ,l, 2 . 

If |A — U 2 I >> equation fll 20 p implies that ‘Fj(A) = 



( 120 ) 


- 2^iti(A-V2-ie) Ir (^2 + *ho(u2 + ^), -Ush + + ^), 

—/i + + <F2('i^2 + ^)j 5 j = 0,l,2, 

( 121 ) 

while, for X — V 2 = y/ti, y = 0(1), then <Fj(A) = 0(1) and 

^j{v2 + y/h) ~ ^ + *ho(u2 + ^), -vsy' + $1(^2 + ^), 

—y' + ^ + $2(^2 + , j = 0 , 1 , 2 . 

( 122 ) 
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Equation fll2ip for |A| >> 1 and fl7S]) finally imply the following result. 

In the space-time regions flll7l) . flll9p . the longtime behavior of the solution 
(/, M, v) of ([2]) reads 


(123) 

^ = -F{-.,t’V-¥i) + o(l), 

where 

-fi + r] + A2X;k,C,v)), 


and the helds Aj{^) = ^j{v 2 + fi/ti), j = 0,1,2 are dehned by the nonlinear 
integral equations 


^jXX)C)'n) 27 ri /r j C) h)) Ch" + 

+ + ^2(P;5,C.>;)), i = 0,1,2. 


(125) 


We observe that the longtime behavior of the solution of (Ej) has the same 
structure as the solution of the linearized version of equation ([2]), except 
from the fact that function F is constructed through the solution of a non¬ 
linear system of integral equations. In addition, the solution, concentrated 
asymptotically on the wave front flllhp . does not break. 

The situation drastically changes if we consider the longtime behavior of 
solutions of (in]),(I7]) in the space regions 

t 2 > 1, X = Vit 2 , y = V 2 t 2 +y, Z = V3.t2, y,Vi,V 2 ,V 3 = 0{1). (126) 

Indeed, in this case, in the region (11261) the integral equations (I115p become 

^ /r + *^o(A'), -'P3(A' - Vi/V3)t2 + <I>l(A'), 

— (A'^ — V2)t2 + y — 2/^2 + ‘^*2(^0^ 5 J = 0; I5 2 


and now the decay of the integral is partially contrasted if ^2 > 0 and i:^/v 2 = 
vi/vs] i.e., on the surface 


± 


y-y 


to 


X 

■) 

z 


^2 

y - —t2 = y 

z^ 


(128) 
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of the (a:, y, z) space. On such a surface, let us consider the case ^/v 2 = Vi/v^, 
it is convenient to make the change of variables y' = (A' — y/v 2 )t 2 in (I127p . 
obtaining, for ^2 S> 1: 

(^) ~ 2 A; /r (v^ + <i>o(v^ + g). -tw' + «!(+ i^), 

— 2^/^fj,' + y — 2/f2 + ‘^*2(v^ ^ is 4 2- 

(129) 

If |A — >> 0(^2 ^)) fhs ^'j ci-re 0(^2 ^md equation fll2Up becomes 


~ - 2nit2{X-ie) (\/^ + MV^ + ^), -^^3/i+ 

+ ■^), -2^y + y - 2 ft 2 + + ■^)y j = 0,1,2, 

while, if A — = y/t 2 , y = 0(1), the $'s are 0(1): 

■s>j(v^+ij) ~hSnii^)Ri{y^^+M^-+ 4.-w'+ 

4i(v^+ ij), -2^2lJ + V- 2/(2 + * 2 (^/ 5 ;+ g)), J = 0,1,2. 


(130) 


(131) 


Equation (I130p for |A| >> 1 and fl76l) hnally imply the following result. 

In the space-time regions (I126p . (11281) the longtime behavior of the solution 
(/,M,u) of ([6]),([7]) reads 


/~^F(f,^,//-f,f2-2/f2) + o(t2l), 
U = -^F(f,^,//-f^f2-2/f2) + o(l), 

_ ‘^^2{fz-\-Uxfx) I ^(t) 


Vx = 


l-2t2fy 


(132) 


where F is dehned by 


F{^,C,V) = -^iI^dyRo(^^ + Ao{y,^,C,v),-Cf^ + 
-2^y + y + A2{y;^,C,v)), 


(133) 


and the helds Aj, j = 0,1, 2 are defined by the nonlinear integral equations 
~ fm ('^ ^o(h' j 'Cl Cl h)) + (134) 

Ali(/i';CiCi^)i-2C/^' + ?? + 242(p';CiCih))i J = 0,1,2. 


By considering the case = —ui/us, we obtain the same result as above. 
Equations fll32p show that, even in the longtime regime, localized solutions of 
equations ([6]) and ([7]) break on the wave front surface fll28p in a way similar 
to the dKP breaking [33]. 
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7 Concluding remarks 

In this paper we have applied the formal Inverse Spectral Transform for 
integrable dispersionless PDEs to the so-called Dunajski hierarchy. We con¬ 
centrated, in particular, i) on the system of PDEs characterizing a general 
anti-self-dual conformal structure in neutral signature, ii) on its hrst com¬ 
muting flow, and hi) on some of their basic and novel reductions. We have 
formally solved their Cauchy problem and used it to construct the longtime 
behavior of solutions, showing, in particular, that unlike the case of soliton 
PDEs, different dispersionless PDEs belonging to the same hierarchy of com¬ 
muting flows evolve in time in very different ways, exhibiting either a smooth 
dynamics or a gradient catastrophe at hnite time. 

In a subsequent paper, we plan i) to investigate the analytical aspects of such 
a wave breaking, taking place, in the longtime regime, on the 3-dimensional 
surface fll28p : ii) to construct, following [17], distinguished classes of implicit 
solutions of the Dunajski hierarchy and of its reductions, concentrating also 
on the class of reductions introduced in [Sj. 
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